Photoemission spectroscopy and sum rules in dilute electron-phonon systems 



P. E. Kornilovitch 

Hewlett-Packard Laboratories, MS 1L-12, 1501 Page Mill Road, Palo Alto, California, 94304 

(February 1, 2008) 



(N 
O 
O 
(N 



A family of exact sum rules for the one-polaron spectral function in the low-density limit is derived. 
An algorithm to calculate energy moments of arbitrary order of the spectral function is presented. 
Explicit expressions are given for the first two moments of a model with general electron-phonon 
interaction, and for the first four moments of the Holstein polaron. The sum rules are linked to 
experiments on momentum-resolved photoemission spectroscopy. The bare electronic dispersion 
and the electron-phonon coupling constant can be extracted from the first and second moments of 
spectrum. The sum rules could serve as constraints in analytical and numerical studies of electron- 
phonon models. 
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Statistical analysis of spectra is an important branch of 
photoemission spectroscopy. It has been known that the 
energy moments of the spectral density of atoms can be 
expressed via the atomic quantum numbers [Q|2|. Such 
integral relations, or sum rules, provide direct link be- 
tween experimental data and the internal state of the 
atoms. Similarly, in crystals analogous sum rules yield 
valuable information about the Bloch states. This is par- 
ticularly true in case of two-dimensional systems where 
the in-plane momentum of the electron does not change 
in the photoemission process. It allows the spectrum to 
be recorded for each in-plane momentum individually. 
In recent years, the angle-resolved photoemission spec- 
troscopy became a powerful experimental technique to 
study the electronic structure of quasi-two-dimensional 
compounds such as the high-temperature superconduc- 
tors and layered colossal magnetoresistive oxides ||-[| . 

This paper discusses statistical analysis of spectra from 
systems with strong electron-phonon (el-ph) interaction. 
In such systems, the mobile carriers deform the lattice 
and form polarons which are bound states of the carriers 
and the lattice deformations. The polarons then define 
the low-energy physics of the system, in particular the 
equilibrium chemical potential. In inverse photoemission 
(IPES), an electron with momentum k is injected into the 
system from outside. Initially it feels the bare electronic 
states in the undeformed lattice and only then begins to 
deform the lattice around itself. Thus it should be ex- 
pected that the energy moments of the IPES intensity 
could provide information about the initial stages of the 
deformation process and therefore about the strength of 
the interaction. It turns out that for low density of po- 
larons the moments of IPES intensity can be calculated 
analytically for a wide class of model el-ph interactions. 
It leads to a series of exact sum rules which the spectral 
density, experimental or theoretical, must satisfy. Below 
I derive the sum rules and discuss their possible applica- 
tions. 



Consider a single band, one optical branch Hamilto- 
nian with an el-ph interaction of arbitrary profile 



H = ^2{e k - M)4 c k - V2fl^2g(n- m)c n c„a 

k nm 



(1) 



Here /i is the chemical potential, x m = Cmy^p is 
the dimensionless displacement of an oscillator at m- 
th lattice site, M and f2 are the mass and the fre- 
quency of the oscillators, respectively, d m = d/dx m , 
c n = iV -1 / 2 ^ k e lkn Ck, and N is the number of lattice 
sites. g(n — m) is the dimensionless coupling constant 
proportional to the force with which an electron at site n 
acts on oscillator m. For simplicity, the fermion spin and 
multiple phonon polarizations are not included explicitly 
in Eq. ([!]) . They are not essential for the purposes of the 
paper and could be added if necessary. The canonical 
Holstein model || follows from the above Hamiltonian 
with the choice g(n — m) = gaS nm . 

Within the sudden approximation, direct and inverse 
photoemission intensities are proportional to the fermion 
spectral functions A_(k, ui) and Af(k, u;), respectively 
(hereafter notation h = 1 is used): 

A_(k, W ) = ^e-^|(j|c k K)| 2 J(a; + ^ -£,-), (2) 



A+fow) = YE e ^IO'l4K)| 2 ^ + ^-^')- (3) 

ij 

Here u> is the energy relative to the chemical potential, 
\i) and \j) are exact eigenstates of TL, with energies £i 
and £j, Z = ^2 n e~^ £i is the grand canonical partition 
function, and /3 — (ksT)^ 1 . The validity of the sudden 
approximation as well as the additional electron-photon 



1 



matrix element, the Coulomb effects and other complica- 
tions, are not going to be discussed in this paper, see e.g. 
M. ft follows from Eqs. (||) and (0) that the two spectral 
functions obey the fundamental relation 



A + (k,u) = e fiu A-(k,u) 



(4) 



Thus, A + can be found if A_ is known with sufficient 
accuracy, and visa versa. In terms of the symmetrized 
function A = A + + A- they are expressed as A_ (k, oj) — 
(e^ + l)" 1 A(k, u) and A+(k, uj) = (e~ 0u + l) -1 vl(k, uj). 

Several theoretical results for the polaron spectral 
function are available. The on-site version of the Hol- 
stein model (e k = 0) is exactly solvable, see e.g. Ref. [ p~0| | . 
For the Holstein model with non-zero kinetic energy the 
spectral function was calculated by Alexandrov and Ran- 
ninger in the Lang-Firsov approximation fll| . Later their 
result was generalized to the long-range interaction (Q) by 
Alexandrov and Sricheewin For the more complex 
Jahn- Teller interaction the polaron spectral function was 
given by Perebeinos and Allen The spectral func- 

tions mentioned are approximate and all have the form 
of a sequence of delta-functions. 

The formalism of statistical moments is best in- 
troduced via the time-dependent correlation functions 
K±(k,t) that are Fourier transforms of A±(k,uj): 



K+(k,t) 



2^ 



e-^ + (k, W ) 



c k e 



-iUtJ \ 



(5) 



/CO J 
^e-^A^cj) 
-co 2n 

= ^J2e-^{j\cl^We- im \j) 



(cle im c k e- mt ) 



(6) 



Taking the m-th derivative with respect to time t and 
setting t — one obtains 



M+(k) 



dw 
2^ 



u m A + {k,uj) 



[ck,«],W]...,«] 4), 



(7) 



m times 



M-(k) 



djj 
2^ 



uj m A(k,Lo) 



J —OO 

(4[[[ck,W] ,H] . . . ,H] 



(8) 



m times 



The moments of the symmetrized function A(k, uj) are 
defined as M m {k) = M+(k) + M m (k). 



The simplest sum rule is the normalization property of 
A(k, uj) which is common to any one-particle fermionic 
spectral function 



M (k) = (c k cl + c k c k ) = 1. 
The first moment of A(k, uj) is given by 

[c k ,7i] = (e k - /^)c k 
V2fl 



(9) 
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e ^ k '-^ n g{n~m)x n c k ,. (10) 



nmk' 



M^k) = ([c k ,H]4 +cl[c k ,H]) 

V2n 



= (£k - M) - 



N 



^ 5 (n-m)(a 



(11) 



The mean displacement of the oscillators (x m ) depends 
only on the total number of carriers in the systems (n), 
and can be calculated as follows. Write the displacement 
operator as x m = (x m ) + y m and substitute in H. Then 
minimization of the free energy F = — \n(e~ l3 ' H ) with 
respect to (x m ) yields 

(x m ) = -(y m ) + V2y^g(n - m)(c] 1 c n ) 

II 

= ^(n)^ fl (n), (12) 



because (y m ) = by definition, and the mean density 
( c n c n) = (n) does not depend on coordinate n. Substi- 
tution in Eq. ( [ll"|) results in: 



Mi(k) = ( £k - M )-2Q(n) 



Es(n) 



(13) 



This is an exact relation for the one-particle spectral 
function. It is a generalization of the sum rule derived 
by Mahan for the on-site polaron jlOj , and by Perebeinos 
and Allen for the Jahn- Teller polaron to non-zero 
polaron density. The main feature of these results is that 
the mean value of the spectral function is not affected by 
interaction with phonons. This conclusion is non-trivial 
in view of significant changes of the spectral function it- 
self, especially in the strong coupling regime. As follows 
from Eq. (13), in the many-body case the first moment 
A^i(k) acquires a correction that is linear in polaron den- 
sity. Proceeding with the calculation one obtains for the 
second moment 

[[<*,«],«] = (4-mH 

V2Q 



N 



e^ k '-^ n g(n-m) 



k'nm 

x [(ek + £k' - 2^).x m + fi<9 m ]c k / 
+ 2 ~W E e 4(k '- k)n 9(n-m) 

k' iimiri 7 

x g(n 



(14) 
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7W 2 (k) = ([[c k ,H],H]cl + ctii^HlH}) = (e k -/i) 2 



N 
2fi 



y^g(n - m)[2(e k - ^)(a; m ) + fi(<9 n 



— 5(n-m)ff(n-m / )(x m a; 



(15) 



Unfortunately, the mean values (x m x m i) cannot be 
calculated analytically which limits application of the 
last expression. In principle, such averages could be 
computed numerically by various techniques, rendering 
Eq. a true sum rule. (It is worth noting in this re- 
spect that in the case of the Hubbard model, the first and 
second moments of the spectral function were calculated 
analytically by Nolting |l4|] . This finding enabled him to 
develop the moment approach to the Hubbard model.) 

Further progress can be made by going to the low 
density limit (n) <C 1. It follows from definitions (|^) 
and (|) that A- on (n) while A + oc 1 — (n). Therefore 
<C M.'tn so that M~ n could be dropped from Eqs. (j|), 
(pd|), and (filf). More importantly, all the phonon oper- 
ators can now be averaged easily. Indeed, in the low- 
density limit the lattice is largely undeformed. The func- 
tion describes the process when an electron is created 
in an empty lattice and removed later. Thus whatever 
phonon operators appear as a result of multiple commu- 
tation with the Hamiltonian, they have to be averaged 
over the system of undeformed independent harmonic 
oscillators, which is always possible. For instance, the 
mean values entering Eq. jl5| ) are (x m ) = (d m ) = 0, 
(x m x m >) = \ coth(i/?ri)(5 mm '. Collecting all the results 
one obtains the following sum rules for the spectral func- 
tion A + (k, us): 



M+(k) = l + 0«n)), 
Mt(k) = (e k - fJ ,) + 0({n)), 
_M+(k) = ( £k -/,) 2 

+ n 2 coth(i/^)]T ff 2 (n) + 0((n)). 



(16) 
(17) 

(18) 



The above exact relations may be used for statistical 
analysis of inverse photoemission spectra (IPES) from 
low-density polaron systems. First of all, Eq. ([l6]) al- 
lows proper normalization of the spectrum for each elec- 
tron momentum k. According to Eq. (|l7|), the first mo- 
ment yields the bare electronic spectrum inside the crys- 
tal. However, in an el-ph system the chemical potential /i 
is lowered by the polaron binding energy E p . As a func- 
tion of k, A4^(k) never gets smaller than E p . In other 
words, one should expect a gap of size E p in IPES. The 
second moment A'l^k) also contains information about 
the strength of el-ph interaction. It is convenient to pa- 
rameterize the latter in terms of the zero order (in hop- 
ping) binding energy of the polaron, = ^X)n5 2 ( n )> 
see [Q. (Ep ^ is employed here as a shorthand nota- 
tion only. It does not include the effects of the kinetic 



energy and therefore it is not the true polaron energy, 
-Ep '' 7^ E p .) The mean squared deviation of the spec- 
trum, (AZTp = Mf - {Mt) 2 = E p Q) ncoth(±f3n), is a 
direct measure of E^ . It has to scale approximately lin- 
early with the IPES gap. Notice, that (Aw) 2 = E p 0) n at 

low temperatures and (Auj) 2 w 2E P °^T at high tempera- 
tures T. 

One should mention that sum rules on the function 
^4_(k,o;) have also been discussed in relation to di- 
rect photoemission. In particular, Randeria et al intro- 
duced an approximate sum rule for the zeroth moment, 
A4q = ?^(k), where n(k) is the mean occupation of the 
state k This relation may not be generalizable to 

polaronic systems because in general n(k) depends on 
the strength of the el-ph interaction and therefore is not 
known a priori. 

Calculation of moments can be continued. The low 
density ensures that any fermionic operator that will arise 
after the multiple commutation of c k with the Hamil- 
tonian, will factor out. The remaining phonon opera- 
tors may be complex but they will have to be averaged 
over the system of free phonons. Such an averaging can 
always be performed. The conclusion is that any mo- 
ment of j4 + (k,o;) can in principle be calculated analyt- 
ically. However, the algebra quickly becomes cumber- 
some. Therefore, I present only two more moments for 
the Holstein model with nearest neighbor hopping, i.e. 
for g(n — m) = goS n m and e k = — J J2i e ' kl ; 1 numbering 



the nearest neighbors: 
M+(k) = (e k -fi) 3 + 



Mj(k) = (e k -^) 4 
+ 2g 2 n 3 (e k 



-2g$Q. 2 {e k -n)+gfa 3 + 0({n)), 

(19) 

-<^ft 2 [3(e k - m) 2 + zJ 2 ] 
[i) + [g 2 +3g 4 T ]n 4 + 0((n)), (20) 



where g\ = g 2 . coth(i/?f2), and z is the number of nearest 
neighbors in the lattice. In deriving (k) and M.~\(\t) 
the following properties of the harmonic oscillator have 



been used: 



i<9m) = -s, and 



= f coth 2 (±/?ft) 



It is interesting that in the fourth moment the hopping 
integral J appears not only via the bare spectrum e k 
but also by itself. Thus A^J(k) is the first moment that 
explicitly distinguishes the dimensionality of the lattice 
through the number of nearest neighbors z. Apparently 
the higher order moments also distinguish different lat- 
tice topologies. 

Apart from their relation to photoemission spec- 
troscopy the sum rules may find more "theoretical" ap- 
plications. First of all, they are valid for the single po- 
laron which continues to draw theoretical interest. For- 
mal transition in the above expressions to the canonical 
ensemble with one particle is done by setting fi — 0. Thus 
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the sum rules could be useful checks for polaron spec- 
tral functions computed numerically by exact diagonal- 
ization (Tcfl , density matrix renormalization group [p/7| , 
quantum Monte Carlo jl8), or any other method. Sec- 
ondly, the sum rules may provide information on the dy- 
namics of polaron formation. For instance, the quantity 
P+(k, t) = |Jf_|_(k, t)\ 2 is the probability that the phonon 
system remains unchanged after time t under the distur- 
bance of an electron with momentum k. On the other 
hand, P+(k,t) = 1 - [MtQc) - Mf 2 (k)}t 2 + o{t 2 ) = 
1 - E^ 0) ncoth(^pn)t 2 + o(t 2 ). It enables one to inter- 



pret r = [Ep 'Q coth(i/3^)] -1 / 2 as the average time of 
phonon emission. This is a non-perturbative result that 
is valid for any model parameters. Finally, the sum rules 
may be employed to improve approximate analytic ex- 
pressions for polaron spectral functions. Consider the 
two functions (fx = 0, T = 0): 



00 2m 

AV ow (k, W ) = (27r) e ^^^-x 
^— ' ml 

m=0 

S(uj + g^Cl — e~ 9 °£ k - mO), 



oo 2m 

^ h (k,c) = (2 7 r)e-o^^_ x 

m— 



(21) 



(22) 



The first function is the one-polaron limit of Alexan- 
drov and Ranninger's result for the Holstcin model in the 
Lang-Firsov approximation jyj. In Eq.(|2l|), g 2 il is the 
zero-order binding energy Ep . A 1 ™ correctly describes 
the low-energy physics of the polaron, most notably the 
polaron band narrowing. However, it satisfies neither of 
the sum rules except the zeroth one. In particular, the 

2 

first moment is equal to e _s °£k instead of the correct 
£k- The second function A!" gh is a "wild guess". It does 
satisfy the first and second moments (but not the third 
and fourth ones) but does not describe the band narrow- 
ing. A_^ g 1 is wrong at low energies but may be accu- 
rate at high energies. This conjecture is supported by 
the anticipation that at high energies the electron should 
become "free" from phonons. Thus a good spectral func- 
tion could be the one that interpolates from A 1 ™ at low 
energies to A^ 11 at high energies. With reasonable pa- 
rameterization, the interpolation parameters could also 
be deduced from the sum rules. 

In conclusion, a series of sum rules on the polaron spec- 
tral function have been derived. The sum rules are valid 
in the low-density limit. The first moment has been gen- 
eralized to finite density. A connection with the angle- 
resolved inverse photoemission spectroscopy has been es- 
tablished. The spectrum should display a gap which size 
correlates with the mean square deviation of the spectral 
function. Several other applications have been identified. 



The sum rules are open to further generalizations to dis- 
persive phonons and more complex interactions. 

I am grateful to A. M. Bratkovsky, D. Stewart, 
J. E. Hirsch, and especially to A. S. Alexandrov for very 
useful discussions on the subject and critical reading of 
the manuscript. 
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